Effects of long range electronic interactions on a one-dimensional electron system 
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The effects of a long range electronic potential on a one dimensional chain of spinless fermions are 
investigated by numerical techniques (Exact Diagonalisation of rings with up to 30 sites comple- 
mented by finite size analysis) and analytic calculations. Due to a competition between the 2kF 
oscillations in the density and the (very slow) log divergence of the long wavelength part of the 
Coulomb potential, the metallic character of the system is enhanced at intermediate (and up to 
quite large) lengthscales. Despite some similarities, we found that this quasi-metallic regime is not 
of the Luttinger Liquid type, as evidenced by strong deviations from predictions of conformal field 
theory and in agreement with the picture of a (very weakly pinned) Wigner crystal. When the 
strength of the Coulomb interaction is substancially increased, we observe a smooth cross-over to a 
strongly localized charge density wave. 

PACS numbers: 71.10.Fd, 71.27.+a, 71.30.+h 



I. INTRODUCTION 

In one dimensional systems long range (LR) Coulomb 
potential plays a particularly important role in electronic 
models [Q . Indeed, as the dimensionality of an electronic 
system is reduced, charge screening become less effec- 
tive in reducing the range of the electron-electron inter- 
actions. The importance of LR interactions between car- 
riers has been highlighted in some experiments on GaAs 
quantum wires or quasi-one-dimensional organic con- 
ductors The role of a 1/r Coulomb repulsion on 
the long distance properties of electrons confined to a 
chain has been theoretically investigated using bosonisa- 
tion techniques ||] . It was found that the charge correla- 
tions decay very slowly with distance suggesting that the 
ground state (GS) is similar to a classical Wigner crystal 
(WC) |^,||. In this analysis quantum fluctuations only 
lead to small oscillations of the particles away from the 
configuration of minimal potential energy. However, this 
previous study was made in the continuum limit and the 
role of the lattice in this context is unclear. Theoretical, 
both analytic ||7|-p^ or numerical [|ri|, and experimental 
results suggest that umklapp processes have a crucial 
importance in ID for the Mott transition. In addition, 
the g — > singular forward scattering interaction and the 
2kp backward scattering might lead to competing effects. 
Although it is believed that, strictly speaking, for a long- 
range 1/r potential of arbitrary magnitude, the lattice 
would always produce a pinning of the WC, such an ef- 
fect is (a) only achieved for very large systems (for not 
too large interaction) due to its logarithmic nature and 
(b) is very weak. Hence, easy de-pinning of the WC can 



occur for mesoscopic systems at small but finite temper- 
ature. 

In order to investigate the interplay between Coulomb 
interaction and umklapp scattering as well as their size 
dependent effects, thermodynamic observables are cal- 
culated on chains with a variety of electron-electron 
interactions using numerical techniques and compared 
to some analytic predictions based on the bosonisation 
method. As the range of the electronic potential in- 
creases, our computational analysis suggests that, at in- 
termediate length scales, the metallic character of the sys- 
tem is enhanced. However, strong umklapp occurs at 
large Coulomb potential leading to an insulating Charge 
Density Wave (CDW) state. Based on our numerical 
data and on analytic calculations we argue in favor of 
a smooth cross-over between a weakly pinned WC (ex- 
hibiting metallic character at length scales and energy 
scales we are dealing with here) and a CDW insulator. 
We believe that the physical behaviors observed here for 
systems of sizes L ^ 20-30 are of particular relevance for 
sub-nanometric quantum wires or carbon nanotubes. 

This article is organised as follows: In Sec. II we inves- 
tigate the role of an extended (but finite) range interac- 
tion. Spectacular oscillations in e.g. the charge stiffness 
are observed numerically as the range of the interaction is 
increased. In Sec. HI, properties of a full 1/r long range 
potential are analysed. First, the finite size behaviors of 
the single particle and two particle charge gaps as well as 
the charge stiffness are studied. Behaviors reminiscent of 
a Luttinger Liquid (with charge gaps oc 1/L and charge 
stiffness oc |1) are seen at small and intermedi- 

ate interactions. Qualitative changes occur at very large 
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Coulomb repulsion with clear signatures of the opening 
of a charge gap and a drastic reduction of the charge stiff- 
ness. Based on a conformal field theory analysis, we then 
analyse the role of umklapp scattering by calculating the 
(size dependent) anomalous exponent. In particular, we 
show that, even at small Coulomb potential, strong devi- 
ations from the Luttinger Liquid (LL) predictions occur. 
Lastly, the momentum dependent low energy spectrum 
as well as the single particle spectral function are com- 
puted confirming the cross-over scenario. 

The model analyzed here consists of a single periodic 
chain of spinless fermions with L-sites and an extended 
interaction. 



H = -t"^{4c^+l + h.c.) 

i 



(1) 



where rii (n) is the local (average) electron density, and 
the rest of the notation is standard (i = 1 is the unit of 
energy). Two types of potentials will be considered: (i) 
a plain LR Coulomb interaction Vi-j — V/\i — j|; (ii) an 
extended range potential, typically Vt-j = V/\i — j\ for 
I* — j| ^ ''max and Vi-j = otherwise, extending up to a 
distance rmax- Using Exact Diagonalization techniques, 
ID closed chains are here analyzed. These numerical re- 
sults are compared to analytic predictions based on the 
bosonisation technique. A uniform ionic background has 
been included in Eq. (||) such that the electrostatic en- 
ergy per unit volume remains finite. In the rest of the 
paper, electron densities n = ^ and n = ^ will be con- 
sidered. Finite size scaling analysis are performed by 
considering rings of typical sizes L — 12, 16, 20, 24 and 
28 [resp. L = 12, 18, 24, and 30] for n = i [resp. n = 5]. 
For an even number of particles closed shell configura- 
tions are used by taking antiperiodic boundary conditions 
(ABC). 



II. FINITE RANGE POTENTIALS 

We start our analysis by first gradually varying the 
range of the electron interaction Vi-j. For a (not too 
large) finite range interaction, the system is expected to 
belong to the LL universallity class characterized by a 
linear collective charge mode (see Appendix). The corre- 
sponding charge velocity can easily be obtained from the 
slope of the low energy mode i.e. the energy difference 
between the first excited state at the smallest finite mo- 
mentum q = ^ and the GS. Numerically, it is found that 
the expected 1 /L^ behavior of the finite size corrections 
I p^ is very well fulfilled so that a simple fit provides ac- 
curate extrapolated results shown in Fig. |l| as a function 
of the potential strength V. It is clear from this pic- 
ture that the leading correction in V is linear and that 
the successive components V2p and V2p+i of the potential 



have quite different effects on the increase of u compared 
to the non-interacting case (u = "Ui? = 2 at half filling). 
Indeed, adding even distance interactions barely affect 
the slopes of u vs V. 




FIG. 1. Extrapolated velocity at density n = 1/2 vs V for 
different ranges rmax = r of the electron-electron repulsion. 
The dotted lines correspond to the theoretical predictions (^) 
for rmax = 1 or 2, 3 or 4, and 5 or 6. 

Some hint for this behavior can be found by using the 
continuum representation explained in Appendix |A|. As 
shown in the Appendix, because of the presence of a 2/ci? 
oscillations in the density, the various range components 
of the interaction interfere destructively. The bare Lut- 
tinger liquid parameters, describing the strength of the 
interactions are affected as (see (|A6|)) 



u/K = VF + -y^ Vr(l - cos(2kFr)) 



(2) 



uK — vp 



For an half filled band 2fcF = tt and vp = 2. In (||) one 
can thus distinguish two contributions: (i) the one com- 
ing from the long wavelength part of the density {q ^ 0) 
for which all interactions simply add. This is the term 
that is responsible for divergence associated with the long 
range Coulomb interaction since '^^Vr diverges logarith- 
mically in this case; (ii) the contribution coming from the 
2kF part of the density. Because the density oscillates, 
interactions of different range can now interfere destruc- 
tively. For long range interactions and large system sizes 
this term is bounded (it corresponds to the Fourier trans- 
form of the interaction potential with a Fourier compo- 
nent of 2kF) and the system is thus totally dominated by 
the long wavelength part. For finite size systems and/or 
finite range interactions the correction coming from the 
2kp term can be sizeable. In particular, at half filling 
one would get 



u^2 + -{2Vi + 2V:i + ---) 

TT 

1-^1(2^1 +21/3 + •••) 
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thus all even terms diseappear in K and u giving a K [u] 
decreasing [increasing] less rapidly than if only the g ~ 
part was concerned. In particular one would have 2V in 
for a nearest neighbor interaction and only ^ 4.36^ 
for a 22 sites interaction. 

The charge stiffness D can be obtained numerically by 
threading a flux $ (in unit of the flux quantum) through 
the ring. The curvature around the minimum of the 
GS energy Eq{^) provides the best estimate, namely 
D — i9^(£'o/i)/90^) [where is a pseudo momen- 
tum (p = i'^'^/L)^]- Note that, for an even number of 
electrons, the minimum occurs at $ — 1/2 i.e. for an- 
tiperiodic boundary conditions corresponding to a non- 
interacting closed shell configuration. Results obtained 
on L = 16, 20, 24 and 28-site rings with Ne ^ L/2 elec- 
trons are shown in Fig. |^. We observe that the charge 
stiffness for Tmax > 1 is systematically enhanced com- 
pared to the case of nearest neighbor (NN) repulsion 
(^max — !)■ Note that for Tmax = 1 a metal- insulator 
transition is known exactly to occur dA, V — 2 so that 
D in Fig. ||(a) [resp. Fig. ||(b)] extrapolates to a fi- 
nite [resp. vanishing] value in that case. An interest- 
ing oscillatory behavior appears as Tmax increases. For 
?'max = 2p -|- 1 (odd), D monotonically increases with 
p. The maximum of the charge stiffness is obtained for 
''max = 2 and D decreases slightly for larger r„iax = 2p' 
(even). The (unique) limit corresponding to a long range 
1/r repulsion would then be reached from below (above) 
when p — > cx) (p' — > oo). We believe that this feature is 
generic and is independent of the magnitude of V as long 
as the potential is short range i.e. r,„ax ^ L. 



D 




0.002 0.004 0.002 0.004 

1/L' 1/L' 



FIG. 2. Charge stiffness at density n = 1/2 vs in- 
verse square length for different ranges rmax = r of the elec- 
tron-electron repulsion with amplitudes (a.) V = 1 and (b) 
V = 4. 

Our numerical results for D can also be explained qual- 
itatively using the bosonisation picture. Indeed, when 
more components Vp of the potential are added, in addi- 
tion to the change of the LL parameter K described by 
( [A^ ) , a dramatic effect occurs also in the strength of the 
umklapp term (see Appendix) 



5 = ^ K cos(2fcFr) (4) 

r 

At half filling, this term is reduced for all even range con- 
tributions (since one has Vi — V2, Vi — V2 + V3 — V4, 
etc. thus enhancing the metallic character of the system 
for intermediate sizes) . Since without such umklapp 
term the stiffness D is not renormalized by interactions 
D — uK — Vp (see (||)), the strength of the umklapp is 
a direct measure of the reduction of the stiffness due to 
interaction effects (the stronger the umklapp term, the 
smaller the stiffness). For a second neighbor interaction 
the umklapp strength g is at its minimum V/2, consis- 
tent with the numerical findings of the largest stiffness. 
Note also that for second neighbor interactions K is not 
decreased compared to the nearest neighbor one, thus 
allowing for a maximum enhancement of the charge stiff- 
ness. For odd range interactions the serie Vi , V1 — V2 + V3, 
14 — ^2 + ^3 — ^4 + 14 and so on is monotonically decreas- 
ing, whereas for even ranges Vi — V2, Vi — V2 + V3 — V4 
and so on is increasing corresponding respectively to in- 
creasing and decreasing charge stiffness in good agree- 
ment with the numerical results. When the whole serie 
is included (for the Coulomb case) the umklapp strength 
is g 0.69V^ and thus much weaker than the nearest 
neighbor one g — V . 

Of course for larger length scales (H) is dominated by 
the logarithmic divergence coming from the g ~ com- 
ponent whereas the 2k p part gives a finite contribution 
and K ~> 0. Because of the presence of the umklapp term 
this always lead to an insulating gapped phase {K < Kc). 
However both for intermediate sizes or finite range inter- 
actions the presence of long range interactions is in fact 
favorable to the metallic behavior. 



III. COULOMB 1/r POTENTIAL 
A. Charge gaps 

We now turn to the case of the full 1/r potential and to 
the investigation of the metal-insulator transition. The 
single particle charge gap at average density fi = is 
defined by Ac,i = Eo{N, + 1) + Eo{N, - 1) - 2Eo{N,). 
Note that the GS energies Eo{Ne ±1) are calculated by 
keeping the ionic charge density constant (i.e. n fixed in 
Hamiltonian (|l|)). Our results for Ac,i at n ~ 1/2 [resp. 
h = 1/3] are shown on Fig. || [resp. Fig. ^. 

At small and intermediate interactions, the charge gap 
seems to vanish as ^ 1/i, a behavior characteristic of 
a Luttinger Liquid (LL). However, we cannot exclude a 
very small extrapolated value. This suggests that the GS 
has a metallic character, at least approximately (i.e. with 
a very small charge gap). On the other hand, at large 
Coulomb repulsion, a charge gap clearly opens up. Note 
that at lower electron commensurability than half filling, 
a significantly larger interaction is needed to clearly see 
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FIG. 3. Single particle charge gap vs 1/L for n = 1/2 and 
several values of the LR Coulomb interaction V. Dotted lines ^^G. 5. Two-particle charge gap vs 1/L for n = 1/2 and 



are guides to the eye. 



several values of the LR Coulomb interaction V . Dotted lines 
are guides to the eye. 



a signature of the gap: V ^ aX n — 1/3 compared to 
only - 5 at n = 1/2. 
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FIG. 4. Single particle charge gap vs 1/L for n = 1/3 and 
several values of the LR Coulomb interaction V . Dotted lines 
are guides to the eye. 



charge excitations (solitons). For example, at half filling 
a dephasing of tt in the 01010101 etc. CDW pattern like 
01011010 etc. produces a charge-1/2 excitation. An ex- 
tra electron introduced in the system will then split into 
two such elementary excitations. Our numerical results 
at large V agree well with the classical estimates e.g. 
Ai^,i = y at n = 1/2 or Ac,i = (1 - at n = 1/3 



We observe that a sizeable gap opens up in the 
LR model only at values of V significantly larger than 
the critical values corresponding to the metal-insulator 
transitions of the related "screened" models where the 
range of the interaction is limited to Tmax — :^ ~ 1 
being the average inter-particle distance). For example, 
at n — 1/2, the single particle gap for a LR potential 
of magnitude V iy > i) is roughly of the same order 
than the gap produced by a short range NN interaction 
of magnitude V /I. This clearly shows again that the ef- 
fect of umklapp scattering is reduced in the presence of 
a LR interaction. However, we believe that, in the case 
of a LR l/r potential, a cross-over rather than a real 
transition occurs for increasing Y . 



In order to estimate the sensitivity of the results to 
shell effects, we have also computed the two-particle 
charge gap at n = 1/2 defined by Ac,2 — Eo{Ne + 2) -I- 
EolNe — 2) — 2i?o(iVe) which only involves even number 
A'e of particles. Results are shown in Fig. ^. For this 
range of parameters and system sizes, we observe, as ex- 
pected, that Ac_2 approximatively equals 2Ac^i although 
the additional finite size correction (typically on 
top of the leading 1 /L term have a different sign. 

Extrapolations of our data for Ac^i to the thermody- 
namic limit have been attempted using a cubic spline 
(including higher powers of 1/L if necessary) and the re- 
sults are shown in Fig. ^. In the classical limit i.e. t = 0, 
electron (hole) excitations are multiple massive fractional 



B. Charge stiffness 

Results for the charge stiffness of the l/r model at 
density n = 1/2 are shown in Fig. |^. Note that we have 
used here two sets of closed rings of sizes L — Ap (12 
to 28 sites) and L = Ap + 2 (14 to 30 sites). ABC and 
PBC are used respectively to ensure {V = 0) closed shell 
configurations in all cases. However, we observe that the 
two sets of data exhibit slightly different scaling behav- 
iors which can be attributed to small shell effects. In 
agreement with the previous results for the single parti- 
cle gap, these data suggest an effective metallic behavior 
at small and intermediate magnitude V while for fairly 
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FIG. 6. Single particle charge gap for n = 1/2 and n = 1/3 
as a function of the magnitude of the LR Coulomb interaction 
V. The dotted lines correspond to the classical behaviors 
Ac,i oc V valid in the t = limit. 



effect is quite small at intermediate coupling V and for 
the system sizes we considered here. Hence, we shall at- 
tempt here a more sensitive analysis based on identities 
from conformal field theories. Note that the latter are a 
priori valid only in the case of standard Luttinger liquids 
(i.e. ID models with short range interactions). There- 
fore, deviations from these expected relations will signal 
non LL behaviors. 

To perform this analysis the charge velocity is needed. 
Results are shown in Fig. The charge velocity increases 
with V and system size as expected. In the continuum 
limit, the WC picture predicts a low energy dispersion of 
the form glog^^^ (q). This logarithmic divergence comes 
directly from the g ~ contribution of the potential to 
u/K as shown in (|^) since Vr is logarithmically diver- 
gent. We thus expect a very weak log^^^ {^/L) divergence 
of u consistent with our data at sufficiently large V. 



large V the charge stiffness is suppressed, even at these 
length scales. Note that, for, let us say, V < A the charge 
stiffness is particularly large, typically of the same order 
as in a non-interacting system. Since D is directly pro- 
portional to the Drude weight in the optical conductivity, 
this signals that most of the optical weight lies at 
(eventhough a tiny pinning frequency might exist). This 
contrasts with the large transfer of weight to finite fre- 
quencies occuring in the insulating phase of the NN t-V 
model for V > 2. 
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FIG. 7. Charge stiffness at density n = 1/2 vs inverse 
square length for a 1/r long range potential of several mag- 
nitudes (as indicated on the picture). Open (closed) symbols 
correspond to an even (odd) number of electrons. 



C. Anomalous exponents 

Although the singular g — > part of the 1 /r potential 
does not alone lead to true localisation, it is nevertheless 
responsible to a if —> LL constant so that umklapp 
scattering becomes relevant and can lead to an insulat- 
ing behavior. However, from the previous analysis, this 




FIG. 8. Charge velocity u in the long range 1/r model at 
density n = 1/2 plotted as a function of the inverse square 
length. Open (closed) symbols correspond to an even (odd) 
number of electrons. 



The anomalous exponent K can be calculated for each 
ring size L by using the conformal identity K = ^. Such 
an id entity can be easily derived from the quadratic form 
( Al) of the boson representation. This quantity is of par- 
ticular importance since, in a conformal invariant theory, 
the charge correlations decay as a power law whose ex- 
ponent is precisely given by K. It will also signal, as 
shown in the appendix, when an umklapp term becomes 
relevant and leads to a gap in the spectrum and an insu- 
lating behavior. 

Our data for K at n — 1/2 obtained for various 
Coulomb repulsion V are plotted vs l/L^ in Fig. ^(a). 
The same set of data is shown in Fig. ||(b) as a function of 
V. For large V, K < Kc so that umklapp scattering is ex- 
pected to be strong in this regime and to produce a large 
charge gap in agreement with the previous findings. The 
results at small and intermediate V are not conclusive 
but suggest that umklapp is weak in this regime. How- 
ever, the (slow) divergence of the charge velocity should 
ultimately lead to K{L) and then, eventually, to the 



5 



opening of a very small charge gap at large sizes. Note 
that K seems to be a more sensitive probe than the sin- 
gle particle gap itself. Indeed, the data in Fig. ||(a) for 
e.g. V = 4: unambiguously show the relevance of umk- 
lapp while the extrapolation of Ac.i in Fig. ||(b) for the 
same parameter is still unconclusive. 



provides strong evidence that the quasi-metallic phase of 
the LR model is not a simple LL fixed point. It would be 
consistent with a system which would be a pinned (by the 
lattice) WC in the infinite size limit and which remains 
metallic of the finite size effects Ifll . 
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FIG. 9. Anomalous exponent K (calculated from D and 
It) of the long range 1/r model at density n = i. (a) Finite 
size behavior for different magnitudes V (as shown on the 
plot). Note that the data obtained on the largest L = 30 ring 
scale slightly differently from the data set obtained on L = 4p 
rings; (b) vs V for different sizes (as indicated in the legend). 
The critical value Kc = 1/2 at half filling is shown by dotted 
lines on each plot. 

One can get additional insights on the possible insta- 
bility of the Luttinger liquid fixed point by considering a 

different estimation K'{L) (see e.g. |15 ) of the anoma- 

^-1 



lous exponent K' — 



where 



is the inverse 



compressibility which can be estimated numerically as 



-\L) 



in order to avoid shell effects. In a 



regular LL the two estimates K and K' converge toward 
the same thermodynamic limit within finite size correc- 
tions ~ . In order to examine the breakdown of the 
LL character, it is then of interest to define the relative 
difference. 



TTU 



(5) 



This quantity is plotted vs V in Figs. |l^(a) and (b) for 
the NN and the long range model respectively using the 
same scale for a direct comparison. As expected, in the 
NN model for y < 2 the relative ratio AK/K is very 
small (typically below 0.01) and strongly increases for 
y > 2 in the insulating CDW phase. The situation is 
drastically different in the LR model where AK/K is 
much larger for all V. Note that the minimum around 
y ~ 6 corresponds in fact to an accidental change of sign 
of the difference K ~ K' and moves continuously towards 
small V for increasing size. We believe that Fig. |l^(b) 




FIG. 10. Relative error ^KjK versus V at density n — ^ 
for a short range NN repulsion (a) and a long range Coulomb 
interaction (b). 



D. Low energy modes and dynamics 

Analysing the structure of the low energy spectrum for 
all q is useful to get further understanding on the quasi- 
metal - insulator cross-over. The dispersion of the GS 
vs momentum is shown in Figs. ^ at constant particle 
density n — 1/m = 1/3 and for two extreme values of 
the Coulomb repulsion, V = 2 and V = 50. Although 
one observes clearly 3 local minima at momenta and 
±27r/3 in both cases, the two spectra are qualitatively 
very different. While at small V linear branches seem to 
appear in the vicinity of the minima, the three almost 
degenerate GS are clearly separated from the rest of the 
spectrum at large V ||l^. This last feature is character- 
istic of translation symmetry breaking occuring in a 2kp 
CDW state in which the (spinless) particles are localised 
at equal spacings. 

A more qualitative analysis can be performed by study- 
ing the size dependence of the gap between the GS and 
the next excited state at momentum q = (so-called 
optical gap) which is shown in Fig. ^(a). At small V, 
Aopt extrapolates for L — > oo to a vanishing or small 
value characteristic of a metallic or quasi-metallic be- 
havior. On the contrary, the extrapolated gap is clearly 
finite for large V as expected in the case of an insulator. 
Note that, for intermediate V , we observe some irregular 
behavior of the data so that no finite size scaling analysis 
was attempted here. Let us point out again that the ab- 
sence of gap for small V is related to the finite size of the 
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FIG. 11. Low energy spectrum of the LR Coulomb model 
at n = 1/3 vs momentum q. At momenta and 27r/3 the first 
excited state is also shown. The definitions of the finite size 
gaps Aopt and A2kp are shown on the plot, (a) V = 2; (b) 
V = 50. 



system and the fact that neither the velocity u nor the 
Luttinger parameter K have been renormalized to zero 
because of the log divergence of 
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FIG. 12. Finite size scaling behaviors of the optical (a) 
and 2kF (b) gaps vs l/L for n — 1/3. 

The instability towards a translation symmetry break- 
ing CDW state can be also evidenced from the energy 
separation A2kp between the absolute GS and the lowest 
energy at momentum 2kF (see Fig. 0(a) for the defini- 
tion). As seen in Fig. [T^(b), at small V /S.2kF scales like 
1/L, a behavior reminiscent of a LL. Expected deviations 
from this behavior due to logarithmic corrections can not 
be seen for such system sizes. For large V, A2fc^ decreases 
faster than 1/L. The data for F = 40 and F = 50 are 
compatible with the expected exponential behavior for 
system sizes larger than the finite correlation length of 
the insulating CDW. 
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FIG. 13. Single particule spectral function vs frequency 
calculated on a L = 20 sites Coulomb ring at half filling. From 
top to bottom, the momenta increase by equal steps from 
to TT. Energies are measured from the chemical potential, (a) 
V = 2; (b) V = 8. 



Lastly, we have computed the single particle spectral 
function A{q,Lo) — — ilmG((7, w), where G{q,Lu) is the 
standard two-points Green function, and results at n = 
1/2 for y = 4 and V = 8 are shown in Figs. |3|(a) and 
(b) respectively. For V = 4 most of the spectral weight 
is located in a single dispersive branch. However, for 
V = 8 when the CDW correlations become stronger, one 
observes clearly a quasi-symmetric "image" branch due 
to Bragg reflection on the CDW (fluctuating) potential. 
These data give a further confirmation of the cross-over 
occuring as the magnitude of the LR Coulomb potential 
increases. 



IV. CONCLUSIONS 

To conclude, we have shown that, for intermediate 
lengthscales an interesting phenomenon occurs for a one 
dimensional system in the presence of long range inter- 
actions. In the infinite size limit, such a system would 
always be insulating due to the log divergence of the long 
wavelength part of the Coulomb potential Vr. How- 
ever, because this divergence is quite slow, for interme- 
diate (and up to quite large !) lengthscales an opposite 
effect occurs due to the 2kF oscillations in the density: 
long range interactions enhance the metallic character 
of the system. Quite generically, umklapp scattering is 
reduced (but not totally suppressed) by Coulomb LR in- 
teractions. At not too large V (typically V" < 20 at 
n = 1/3) the ID electron system exhibits quasi-metallic 
behaviors with vanishing (or tiny) mass gaps in agree- 
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ment with previous findings on spinfuU LR models 
This phase, however, shows deviations from conformal 
field theory predictions. A physical description in term 
of weakly pinned WC is proposed. At larger Coulomb 
interaction a cross-over towards a strongly localised 2kp 
CDW occurs. Analytical calculations based on bosonisa- 
tion techniques confirm our numerical findings. 

We thank IDRIS (Orsay) for allocation of CPU time on 
the C94 and C98 CRAY supercomputers. TG would like 
to thank the ITP where part of this work was completed 
for support (under NSF grant No. PHY94-07194) and 
hospitality. 

Note added: While completing this work w e be- 
cam e aware o f related work by G. Fano et al. |cond- 



ma t l/990914Ctl in agreement with our analysis in the small 
V regime. 
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APPENDIX A: BOSONISATION 

One can study the effect of Coulomb interations on the 
Id electron gas by using a boson representation of the 
fermion operators (bosonisation) [ p^[pj| ]. Introducing a 
field (f) to represent the long wavelength fiuctuations of 
the density, the kinetic energy in (|^) becomes 



If u 
H = — dxuK{TrIlf + — (V0)^ 
zvr / K 



(Al) 



where (f> and 11 are canonically conjugate operators. For 
the free {V = 0) Hamiltonian u — vp and K = 1. The 
parameters u and K are the only one necessary to de- 
scribe the low energy properties of the metallic phase and 
are known as Luttinger liquid parameters. They can be 
computed in perturbation in the interaction V by using 
the expression of the density 



1 



p{x) = po V0 + 



1 



J{2kFX+24>) 



(A2) 



TT 27ra 

where a is a short distance cutoff that can be identified 



with the lattice spacing. Using (A2) for the interaction 
term in (|^) one gets 



i?int=5]^.-,{^(V0.)(V0,) 



(A3) 



(27ra)' 



:2cos(2/ci?(a;i — Xj) + 20^ — 20^)} . 



By considering that the fields (j) are slowly varying at 
the scale of the lattice one can expand (A3). To expand 
safely the fields should be normal ordered using 



gi2(</,(ri)-0(r2)) ^. gi2(0(ri)-0(r2) . i ([2(0(ri )-</.(r2 ))]" 



At the lowest order in V, (A4) becomes 



gi2(0(ri)-</.(r2)) ^. gi2(0(ri)-0(r2) . 



(ri - 7-2)^ 



(A4) 



(A5) 



In the normal product one can expand (/)(ri) — 0(r2) ~ 
(ri —r2)V(j){R) where R = (v i + r2)/2. One thus recovers 
a quadratic form similar to (Al) but with 



2 X - 

u/K = VF + -} Vr(l - cos(2kFr)) 



(A6) 



uK — vp 



where r — 1,2,3,... is the distance on the lattice. The 
above expression is valid at the lowest order in V . For 
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finite V the various coefficients in it will get renormalized 
by irrelevant operators. For a long range potential the 
g — > part is dominant and gives rise to the logarithmic 
divergence of Aq responsible for the one dimensional plas- 
mon. Due to the 2kp oscillations in density the 2fci? com- 
ponent of the potential also occurs. Note that although 
( [A^ ) looks superficially similar to what one could obtain 
using g-ology arguments, such arguments would naively 
suggest that gi = — 52- This would lead the V{2kp) part 
to enter both the u/K and uK term and lead to an in- 
correct renormalization of the uK term. Such a change 
of uK is impossible for a well defined lattive model [gOj , 
where one should have uK = vp as a consequence of 
Galilean invariance in the absence of an underlying lat- 
tice. 

In the presence of a lattice one should take into account 
that the momentum is only conserved modulo a vector 
of the reciprocal lattice (i.e. for example for half fi lling 



Akp = 27r). When taking this effect into account (A3) 
gives also rise to the so-called umklapp terms. These 
terms are the one responsible for the Mott transition and 
the opening of the charge gap ||] . They can also be com- 
puted, and read e.g. for the half filled case 

-ffumk = TTT^ / da;cos(4(/)(j;)) (A7) 
{2T:aY J 

where 

g = ^VrCos{2kFr) . (A8) 



Taken together with ( |Al[ ), ( |A7] ) is able to open a gap 
in the spectrum and lead to a Mott insulator (or equiv- 
alently to a Wigner crystal pinned on the lattice) for 
K < Kc, where Kc depends on the commensurability. 
For half fiUing and would be = 1/(2^2) for 

a commensurability of order n. 
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